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ABSTRACT 

We propose new diagnostics for circumstellar interaction in Type IIP super- 
novae (SNe IIP) by the detection of high velocity (HV) absorption features in 
Ha and He I 10830 A lines during the photospheric stage. To demonstrate the 
method, we compute the ionization and excitation of H and He in supernova 
ejecta taking into account time-dependent effects and X-ray irradiation. We find 
that the interaction with a typical red supergiant wind should result in the en- 
hanced excitation of the outer layers of unshocked ejecta and the emergence of 
corresponding HV absorption, i.e. a depression in the blue absorption wing of Ha 
and a pronounced absorption of He I 10830 A at a radial velocity of about —10^ 
km s~^. We identify HV absorption in Ha and He I 10830 A lines of SN 1999em 
and in Ha of SN 2004dj as being due to this effect. The derived mass loss rate 
is close to 10~^ Mq yr~^ for both supernovae, assuming a wind velocity 10 km 
s~^. We argue that, in addition to the HV absorption formed in the unshocked 
ejecta, spectra of SN 2004dj and SN 1999em show a HV notch feature that is 
formed in the cool dense shell (CDS) modified by the Rayleigh- Taylor instability. 
The CDS results from both shock breakout and radiative cooling of gas that has 
passed through the reverse shock wave. The notch becomes dominant in the 
HV absorption during the late photospheric phase, > 60 d. The wind density 
deduced from the velocity of the CDS is consistent with the wind density found 
from the HV absorption produced by unshocked ejecta. 

Subject headings: stars: mass loss — supernovae: general — supernovae: indi- 
vidual (SN 1999em, SN 2004dj) 
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INTRODUCTION 



The link between a certain Type IIP supernova (SN IIP) and its main-sequence pro- 
genitor is poorly known despite the widely accepted view th at these supernovae arise from 



the initial mass range of 10 — 25 Mq (e.g., iHeger et al.ll2003l ). The primary reasons for the 



uncertainty in this field are the small number of presupernova (pre-SN) mass determinations 
and the uncertain amount of matter lost through winds, presumably dominated by the red 
supergiant (RSG) wind. While the former is becoming well-studied, the latter uncertainty 
is related to our poor knowledge of the complicated physics of mass loss and the unsatis- 
factory situation with the empirical measurements of mass loss from RSGs. Even for the 
most studied close massive RS G, Betelgeuse, the range of observational estim ates is large: 
2 X 10"^ - 1.5 X 10"^ Mq yr"^ J Jura fc KleinmannI Eoool : Ivan Loon etal1l2005h . 



The situation for SNe IIP has some promise because the dominant mass loss at the RSG 
stage may be observed by the detection of circumstellar (CS) interaction at X- ray and radio 



wave lengths. The Type IIP SN 1999em was detected in X-rays with Chandra (jPooley et al. 
20021 ). leading to the mass loss estimate M (1 — 2) x 10~^ Mq yr~^ (assuming a wind 
velocity of 10 km s^^). A recent study of available X-ray and radio data for S Ne IIP produced 
a ran ge of mass loss rates of pre-SN IIP M ~ (1 — 10) x 10^^ Mq yr~^ ( iChevalier et al. 



20061 ). The application of these rates to the full RSG stage (~ 10^ yr) suggests the loss 
of 1 — 10 Mq for pre-SNe IIP. The wide range of estimates emphasizes the need for the 
individual determination of the wind density for each particular SN IIP under consideration. 
Unfortunately, this is not always possible since X-ray and radio observations of SNe IIP are 
often not available. 

Here, we propose two new diagnostics for the wind density in SNe IIP that could help. 
Both rely on spectroscopic observations of Ha and He I 10830 A at the photospheric epoch. 
The first one is based on the fact that the interaction of SN ejecta with the wind results 
in the emission of X-rays from both forward and reverse shocks. The X-rays cause ioniza- 
tion and excitation of SN ej ecta that may be revealed, e.g., through specific emission lines 
( 1 Chevalier fc Franssonlll994l ). Unfortunately, in SNe IIP the wind density is low and emis- 
sion lines caused by CS interaction are extremely weak and cannot be detected. We find, 
however, that the excitation of H and He produced by X-rays in SN IIP ejecta turns out to 
be sufficient to be detected as high velocity (HV) absorption features in Ha and He I 10830 
A lines against the bright SN IIP photosphere. This is the core of our proposed diagnostic 
for the wind in SNe IIP. The second proposed probe for CS interaction exploits the possi- 
bility that a cool dense shell (CDS) might form at the SN/CS interface because of radiative 
cooling. The CDS excited by X-rays could become visible as narrow HV absorption in Ha. 
The velocity of this absorption would provide a direct measure of the expansion velocity of 
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the SN/CS interface, a valuable dynamical characteristic of the CS interaction. 



The identification of th e expected HV lines iii observed spectra is complicated by the 
presence of weak metal lines (jPessart fc Hillierll2005l ) . However, we argue that HV lines of Ha 



and H;^ have been observed in spectra of SNe HP, as previously discussed by [Leonard et al. 
J2002ah for SN 1999em. 



The paper is organized as follows. We start with a description of the interaction model 
and CS interaction effects in the Ha and He I 10830 A absorption lines formed in the un- 
shocked ejecta during the photospheric epoch (§[2]). We then compare our CS interaction 
models with the available spectra of SNe HP in the Ha and He I 10830 A lines and estimate 
the wind density for particular SNe HP (§[3]). In §|3]we address the issue of the Ha absorp- 
tion in the cool dense shell at the SN/CS interface of SNe HP. We discuss implications of 
our models in the last section. 



2. MODEL 

The model for ejecta-wind interaction effects in Ha and He I 10830 A lines formed 
in the SN ejecta consists of three major parts: (i) an interaction model that provides the 
dynamical evolution of the SN/wind interface and the X-ray emission from the reverse and 
forward shocks; (ii) a model for the ionization and excitation of H and He in the unshocked 
SN ejecta irradiated by X-rays; and (iii) the calculation of line profiles. We perform (iii) 
using either a standard Sobolev approximation or the direct integration of the radiation 
transfer equation, depending on the validity of the Sobolev approximation. 



2.1. The interaction model 



The interact i on of SN ejecta w ith CS wind leads to a canonical double-shock structure 
(jChevalierl Il982al : iNadyozhiru Il985l ) with the forward shock propagating in the CS gas and 
the reverse sho ck in the SN ejec ta. We treat the CS interaction of ejecta in the thin shell 
approximation (jChevalieiill982bl ) in which the double-shock layer is reduced to an infinitely 
thin shell. We assume that the freely expanding {v = r/t) SN envelope has a sharp boundary 
at the velocity Vc and begins to interact with a wind starting some moment (t ~ 1 day) which 
corresponds roughly to the shock breakout phase. Free expansion is expected to take several 
doubling times of the initial radius, or several days, to be set up, but the observations we 
are modeling are at later times. The maximum velocity Vc is set by the escape of radiation 
from the shock wave at shock breakout. 
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For the density distribution p{v) in the SN IIP envelope we use an analytical expression, 



Po 



1 + [v/voY 



that closely approximates the c ombination of an inner plateau and outer power law tail found 
in hydrodynamic models (e.g.. lUtrobinI 120071 ) . The parameters po and Vq are determined by 
the kinetic energy ejecta mass M, and k: 
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9. We 



adopt the boundary velocity for a typical SN IIP of Vn = 1 . 5 x 10^ km s ^, which is consistent 



with models of shock breakout (IMatzner fc McKed Il999l : IChevalier et al.ll2006l ). Also, the 
blue edge of the H/3 e mission in SN 1999gi on day 1 was observed to be at 15,000 km s~^ 
Jleonard et al.ll2002bh . 



Hydrodynamical modeling predicts that at the shock breakout phase a thin dense shell 



forms at the outer boundary (IGrasberg et al.l Il97ll : iFalk fc ArnettI 119771 : IChevalierl Il981 



UtrobinI 120071 ) . The physics of the 'boundary shell' formation is in the transition from the 
adiabatic to the radiative regime of shock wave propagati on in the outerm ost layers of the 
exploding star. Simple considerations of radiative diffusion (j Che valierl 1 1 98 if ) give an estimate 
for the shell mass: 



M, ^ 2 X 10 
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0.4 cm^ g 
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(4) 



where R is the progenitor star radius and Vb is the velocity at shock breakout, before free 
expansion has been established. During the acceleration phas e the boundary shell is subject 
to the Rayleigh- Taylor (RT) instability (IFalk fc ArnettI 119771 ) and therefore could be corru- 
gated or even fragmented. However, it is not clear whether the RT instability results in full 
fragmentation; 2D or 3D hydrodynamic modeling of the phenomenon is lacking. We will 
treat the boundary shell in our model as an intact thin spherical shell which is already in 
place at the initial moment. The mass of the boundary shell is presumably equal to the 
mass of the outer envelope {v > v^) with an extrapolated power law p oc . This provides 
a good estimate of th e shell niass fo r SN 1999em, ~ 3 x 10"'^ M©, in agreement with hydro- 
dynamical modeling (jUtrobinl 120071 ). The boundary shell is considered below as a seed for 
the CDS that may, or may not, grow further due to radiative cooling at the reverse shock. 
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To calculate the X-ray emission from the reverse and forward shocks we assume that the 
postshock layers are uniform and their densities are 4 times larger than the corresponding 
preshock density. We allow the electron temperature to be partially equilibrated, assuming 
the following interpolation: 



max 



t 



-T 

J. p, 



(5) 



where teq is the electron-ion temperature equilibration time (ISpitzerl Il962l ) and Teq is the 



equilibrium shock temperature with Tg = Tj. The lower limit Tg/Teq = 0.1 is in dicated by 
obser vational data on supernova remnants with shock velocities of ~ 1000 km s~^ (IRakowski 
20051 ) . The reverse shock is equilibrated during the photospheric phase for a mass loss rate 
of ~ 10~^ Mq yr~^ with a wind velocity of 10 km s~^, while the forward shock is not. Since 
the reverse shock is equilibrated and dominates the X-ray luminosity, the uncertainty in Tg 
does not affect our results. The shock X-ray luminosity is expressed in terms of the kinetic 
luminosity as Lx = r/x-^kin- The efficiency parameter r/x is smoothly interpolated between 
extreme cases, r/x = t/it + t n), where tr. is the coo l ing ti me and is calculated for each shock 
using the cooling function of [Sutherland fc Dopital (Il993l ). The wind is assumed to be steady 
with the density p oc r~^. Since all the interaction effects for low wind velocities [u) are 
determined by the wind density p oc M/u, we use hereafter the dimensionless wind density 
parameter w = M-g/wio, where M_6 is the mass loss rate in units of 10~^ Mq yr~^ and uio 
is the wind velocity in units of 10 km s~^. 

The interaction model for a SN IIP with the wind density w = 1 is illustrated in Fig. 
[U The two cases shown have the same ejecta mass M = 19 Mq but different kinetic energy: 
E = 1.3 X 10^^ erg and = 0.65 x 10^^ erg. The high energy case corresponds presumably 
to SN 1999em (jUtrobinll2007l ). while the low energy case illustrates SN 2004dj, as we will see 
below. Accordingly, we adopt Vc = 15, 000 km s~^ in the high energy case and Vc = 13, 000 
km s~^ in the low energy case. The model with parameters for SN 1999em reproduces 
the observed (absor bed) Chandra X-ray luminosity of SN 1999ern in th e 0.5 — 8 keV band 
(jPooley et al.ll2002l ) for a distance D = 11.7 Mpc ([Leonard et al.ll2003l ) quite satisfactorily 
(Fig. [T]), which indicates that w = 1 is a reasona ble choice for the wind density. The 



estimated Galactic absorption, iVn = 6 x 10 cm (jPooley et al.ll2002l ). is small and does 



not significantly affect the observed X-ray luminosity. The X-ray luminosity from the reverse 
shock dominates in both cases. The typical electron temperature for the reverse shock at 
t > 20 d is ~ 1 keV, while for the forward shock it is ~ 10 keV. The early behavior of 



the velocities differs from the self-similar evolution (jChevalierl Il982bl ) because of the ejecta 
density cut-off in our model and the presence of the boundary shell at the initial moment. 
The mass of the boundary shell is 3 x 10^"^ Mq and 1.1 x 10~^ Mq in the high and low energy 
models, respectively. The total mass of the CDS attained at the time of the termination 
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of the radiative regime in the reverse shock (trad) is 6 X 10"^ Mq and 3.2 x 10"'' Mq with 
^rad = 21 d and trad = 18 d for the high and low energy models, respectively. Remarkably, 
the reverse shock is adiabatic in the self-similar models for the same parameters through this 
epoch. This fact is related to the different dynamics of the self-similar model at an early 
stage, which results in a lower preshock density at the reverse shock. 

The spectrum of X-rays is described below using a standard approximation Fe oc 
exp(— E'/Tg). We ado pt q = 0. 4 for the forward shock, a reasonable approximation 



for temperature ~ 10 keV (jCox 2000), while for the low temperature reverse shock we take 
q = I- The latter choice qualitatively describes the calculated energy distribution of X-rays 
with a significant co ntribution of emission lines for the electron temperature Tg = 1 keV 
jNvmark et allbood ). 



2.2. Model of H and He excitation 



The hydrogen and helium ionization in a SN IIP a tmosphere at the photospheric epoch 
is time-dependent (lUtrobin fc Chugaill2002l . |2003| . l2005l ). To take this effect into account, we 
use a hydrogen atom represented by four levels and continuum. All the essential radiative 
and coUisional processes for hydrogen are taken into account, including non-thermal ioniza- 
tion and excitation induced by X-ray absorption. The excitation of the helium 2^S level 
is determined by a simple model with two excited levels (2^S and 2^S). We consider only 
major processes for the population of the 2^S level: nonthermal excitation and ionization 
with subsequent recombination into triplet st ates, depopula tion by the collisional deexci- 
tation transition 2^S-2^S, Penning ionization (jChugail Il99ll ) and continuum absorption in 
the 10830 A line (2^S-2^P transition) with the subsequent spontaneous transition 2^P-1^S 
terminated by photon escape. The energy balance includes the energy deposition due to the 
absorption of X-rays and the hydrogen photoionization from excited levels by photospheric 
radiation. The cooling includes adiabatic loss {pdV work), free-free, free-bound, bound- 
bound emission of hydrogen, and cooling in the [01] 6300, 6364 A and Mgll 2800 A lines. 
The Fe II line cooling is approximately included by multiplying the Mgll 2800 A cooling 
rate by a factor of two. Given the possibility of a high temperature regime (Te > 20, 000 K) 
in the oute r layers, the cooling at this range of temperatures is treated using the cooling 
function of [Sutherland fc Dopital (119931 ). 



The distribution of the ionization and excitation of H and He I in the atmosphere at one 
time is calculated by tracking the evolution of a Lagrangian zone after its emergence from the 
photosphere. The initial ionization and excitation at the photosphere is set assuming Saha- 
Boltzmann equations for the effective temperature. Subsequent ionization and excitation are 
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computed by solving the system of time-dependent kinetic equations together with the time- 
dependent energy balance. The average intensity of the photospheric radiation is assumed 
to be Ji, = Wly, where W is the dilution factor and 1^ is the photospheric brightness. The 
luminosity evolution at the photospheric epoch {t < 100 d) is described by the expression 



L = Lo[l + exp(-t/ti)]. 



(6) 



With Lq = 10^2 erg s ^ and U = 10 d, this corresponds to the bolometric luminosity of 



SN 1999em (lElmhamdi et al.ll2003l ). assuming a distance of 11.7 Mpc. The photospheric 
radius is determined by the veloci ty at the photosphere and the age (rp = Vpt), where Vp is 



deduced from observational data of lLeonard et al.l (l2002al ). The spectrum of the photospheric 
continuum is taken to be blackbody with T = Tgfj in the wavelength region longward of the 
Balmer jump, while the ultraviolet (UV) radiation shortward of the Balmer continuum edge is 
suppressed by a factor /uv < 1 which was assumed to monotonically drop after the explosion 
from unity to 0.03 on a timescale of 1 d. This evolution was deduced from the behavior 
observed in SN 1987A (jPun et al.lll995l ). We found that the results are only weakly sensitive 
to the specific form of the variation of fyy. The energy deposition rate in the atmosphere 
due to the absorption of X-rays is calculated with an absorption coefficient kx = lOOE^^^'^ 
cm^ g^i (where E is in keV). The deposi ted energy is shared betw een h eating, ionizatio n, 



and excitation according to the recipes of iKozma &: FranssonI (119921 ) and IXu et al.l (Il992l ). 



The solution of the kinetic equations for each Lagrangian zone provides populations of 
the hydrogen levels n2 and n^, and the 2'^S level of He I in the atmosphere. These values are 
then used to calculate line profiles of Ha and He I 108 30 A. We compared our m odel in the 
case of zero wind density with the detailed model of lUtrobin fc Chugail (120051 ) and found 
that the present model results in slightly weaker Ha absorption; at the relative intensity of 
0.9, the radial velocity in the simple model is lower by ~ 500 km s~^. This drawback of 
the simple model does not noticeably affect the amplitude of the calculated HV absorption 
formed in the external layers, where the ionization unaffected by X-ray absorption is small. 



3. HIGH VELOCITY FEATURES IN H AND He LINES 

3.1. Model properties 

The model line profiles of Ha and He 1 10830 A are computed for different wind densities, 
ages, ejecta mass and energy (Table 1). The standard model is characterized by M = 19 Mq, 
E = 1.3 X 10^^ erg, k = 9. In all cases Vc = 15, 000 km s~^ with the exception of the very 
low energy model (eW) for which Vc = 6500 km s~^. We designate a model at a particular 
time by adding the age in days to the model name; e.g., W50 stands for model W on day 
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50. The effect of the wind density on Ha on day 50 is shown in Fig. [2]a for the wind density 
parameters w = 0.5, 1, 2, and 4. Without a wind, the Ha absorption on day 50 forms because 
of the time-dependent ionization effect; a high ionization is maintained due to Lya trapping 
and reionization from excited levels of hydrogen. The CS interaction results in the ionization 
and excitation of the outer recombined layers of unshocked ejecta and, as a result, a HV 
absorption feature appears in the blue wing of the undisturbed profile. Fig. [3] illustrates 
how Ha absorption is modified by the CS interaction. For w = 0.5, the CS interaction effect 
is small (Fig. [2]a) and the wind density w = 0.5 seems to be the least that could be detected 
by HV absorption in Ha. For w ^ 1, the strength of the HV absorption increases with the 
wind density roughly as w"^ and for w > 2 becomes saturated. The decrease of the velocity 
of the blue absorption edge is due to the deceleration caused by the dense wind. The fact 
that for w = 4 the HV absorption merges with undisturbed Ha absorption implies that in 
the case of a very dense wind it would be difficult to infer the density parameter w using 
only the Ha absorption. 

The evolution of the HV absorption (Fig. [2]6) shows that there is an optimal phase at 
about the middle of the plateau (~ 50 d) when the CS interaction effect is most pronounced. 
At an early stage (e.g., on day 20) the CS interaction effect is not clearly seen because of 
the merging of HV absorption with the strong undisturbed absorption, while at a late stage 
(e.g., 80 d) HV absorption becomes very faint. The sensitivity of HV absorption to the 
ejecta density is shown by the model mW50 with a mass 10 Mq and energy 9 x 10^° erg. 
The energy is chosen to match the blue wings of the main absorption of models mW50 and 
W50. The HV absorption in the low mass case turns out the same as in the model W50, 
i.e., the effect of reducing the ejecta density by a factor of 2 is negligible. We also consider 
the CS interaction effect for SNe IIP with a very low kinetic energy, E = 3 x 10^° erg. 



which presumably applies to the low- luminosity SNe IIP, e.g., SN 1999br (jPastorello et al. 



20041 ). This model with Vc = 6500 km s~^ and w = 1 predicts strong HV absorption which is 
merged with the undisturbed absorption (Fig. [2](i). For all the models when HV absorption 
is unsaturated, the Sobolev optical depth of this feature for w ~ 1 is r < 1. This implies 
that in the H/3 line the optical depth < 0.14, so we do not expect that HV absorption in this 
line will be pronounced (but see § S]) . 

We examined the sensitivity of HV absorption to the power law index q of the X-ray 
spectrum of the reverse shock. We found that a softer spectrum {q = 1.2) and harder 
spectrum (g = 0.8) produce slightly deeper HV absorption, so the preferred case q = I 
corresponds approximately to a minimal effect of the CS interaction on Ha. We checked the 
sensitivity of HV absorption to the adopted SN bolometric luminosity as well. The HV ab- 
sorption strength anticorrelates with the bolometric luminosity, but the dependence is weak. 
The mechanism for this dependence is the depopulation of excited levels by photoionization. 
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The CS interaction effect in the He I 10830 A hne substantially differs from that in Ha 
in one important respect: in the absence of a wind the model does not predict noticeable 
absorption in He I 10830 A (Fig. |l]a), except for a very early stage {t < 10 d). The wind 
density w = 0.5 produces weak HV absorption, on the verge of detectability, while the case 
w = 1 yields HV absorption with a relative depth of ~ 0.2 and could be easily detected. The 
absorption gets markedly stronger for larger w, showing the sensitivity of He I absorption 
to the wind density in this particular range of the w parameter. Remarkably, while in Ha 
the contribution of the interaction effect for w = 4 is difficult to discern (Fig. [2]a), the He I 
10830 A HV absorption provides a probe for such a dense wind. 

The strength of HV absorption in He I 10830 A increases between days 20 and 50 but 
gets weak on day 80 (Fig. Wi>), resembling the behavior of HV absorption in Ha. The effects 
of the ejecta density (Fig. Hie) are similar to those in Ha. The HV absorption for the low 
energy case is very strong and easily discernable, unlike in Ha. 



3.2. Comparison with observations 



3.2.1. HV absorption in Ha 



Spectra of SN 1999em (lElmhamdi et al. I l2003h on days 41 and 54 (assuming the explosion 
date JD 2451476) show a depression in the blue wing of the Ha line (Fig. [5]) at the predicted 
position of HV absorption (—10,000 to —12,000 km s~^). To model the Ha profiles we 
adopted the standard parameters M = 19 Mq, E = 1.3 x 10^^ erg, k = 9, and Vc = 15, 000 
km s~^. We find the best value of the wind density parameter is w = 1 with an uncertainty 
of ~ ±0.07, assuming the other parameters to be fixed. The indicated error reflects only the 
accuracy of estimating this parameter from the best fit and does not include errors in the 
observed spectrum and uncertainties in the model. Also on day 54, the model reproduces 
the general strength of the observed HV absorption (Fig. Yet we find that the observed 
HV absorption on day 54 seems to contain a notch component that is deeper than in the 
model. 

The question arises of why the model fails to describe in detail the undisturbed Ha. 
The model emissio n is weaker, while the absorption is stronger, than is observed (Fig. [5] 
and lUtrobinI 120071 ). The weak model emission is probably related to our c rude description 
of the UV continuum in the SN atmosphere (cf. lUtrobin fc Chugail l2005l ) . However, this 
cannot explain why the observed absorption is shallow. In fact, a check of available spectra 
of SNe IIP indicates that the Ha absorption seems to be shallow in any normal SN IIP (e.g., 
SN 2004dj, see below). In contrast, for the peculiar SN 1987A related to the explosion of 
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a blue supergiant, the Ha abso rption is deep and well reproduced by the time-dependent 



model (lUtrobin fc Chugaill2005l ). We speculate that this is a strong indication that the shal- 
lowness of the Ha absorption in normal SNe IIP is related to the red supergiant structure 
of the presupernova. For instance, convection in the red supergiant atmosphere could pro- 
duce large amplitude density perturbations that are significantly amplified during the blast 
wave propagation. As a result, the outer layers of ejecta would acquire a clumpy structure 
responsible for the shallowness of the Ha absorption. 

Return ing to the HV absorpti on, in Fig. [6] we show the model fit to the observed Ha in 



SN 2004dj flKorcakova et al.l 120051 ) on days 55 and 64 (for an explosion date JD 2453185). 
The model parameter set is M = 19 M©, A; = 9, as for SN 1999em, but E = 6.7 x 10^° 
erg. One could retain the energy 1.3 x 10^^ erg used for SN 1999em and vary M. However, 
in this case one requires M ^ 50 Mq to attain a comparable fit. Since the ejecta mass 
unlikely exceeds 20 Mq, the lower energy case for SN 2004dj is preferred. In accord with 
the lower energy we take Vc = 13, 000 km s~^ and lower velocities at the photosphere by a 
factor 0.9 compared to SN 1999em. The derived wind density {w = 1.27±0.1) is determined 
by the strength of the HV absorption and its position, while the lower energy is needed 
to reproduce the correct position of the HV absorption. To check the effect of the adopted 
boundary velocity Vc, we modeled the case Vc = 11, 000 km s~^. In this case the wind density 
is lower, w = 0.9. We prefer Vc = 13, 000 km s~^ on the basis of the arguments presented 
below in Sec. 



On day 55 the HV absorption in the observed spectrum has the appearance of a relatively 
broad shallow depression similar to SN 1999em on day 41. A bit later, on day 64, a notch 
("V-shaped" narrow absorption) appears at a radial velocity ~ —8200 km s~^ (Fig. [6]). 
The feature dominates the HV absorption until at least day 102 in spectra of SN 2004dj 



reported by lVinko et al.l (120061 ). Neither the strength nor the shape of the notch at the late 
photospheric stage t > 60 d can be produced by a model in which the HV absorption forms 
only in the unshocked ejecta. This suggests that the notch in Ha has a different origin which 
is addressed in § HI 

The HV absorption in SN 2004dj has a markedly lower velocity compared to SN 1999em 
at a similar epoch (8200 vs. 11,500 km s^^). This fact is significant for the origin of the 
depression in the blue wing of Ha since it rules out the possibility that the feature might be 
produced by some unidentified metal line. Photospheric velocities at this epoch are similar 
in both SNe within several 100 km s~^, so the expected positions of a metal line absorption 
should coincide in both SNe within the same range of velocities. Since this is not the case, 
we conclude that HV absorption in these supernovae cannot he produced by an unidentified 
line. An additional argument against the 'metal line' option is that SN 1987A does not show 
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HV absorption in Ha at all. This would be difficult to understand in the context of the metal 
line option, but is easily expla ined by the very low density of the pre-SN wind (ty ~ 2 x 10^^, 



Chevalier fc Dwarkadaslll995l ) that is incapable of producing detectable interaction effects. 



3. 2.2. HV absorption m Hel 1 0830 A 



The HV absorption of He 1 10830 A in SN 1999em on day 20 was identified by lElmhamdi et al, 

(120031 ) and attributed originally to radioactive excitation by external ^^Ni. Here we abandon 



this explanation and propose excitation of He in the external lay ers by CS interaction. The 



He I 10830 A line of SN 1999em on day 20 flElmhamdi et al.ll2003[ ) is shown in Fig. [7| together 



with models for w = 1 and the best value w = 1.15, which has an uncertainty of ±0.07. The 
wind density parameter derived from the HV absorption strength in Hel 10830 A agrees 
well with the value w = 1 found from the Ha line. The clear-cut observed profile and the 
sensitivity of HV absorption of He I 10830 A to the wind density makes this line a valuable 
diagnostic tool for probing the CS wind in SNe IIP. 

Interest ingly, based on a time- dependent model for a SN IIP spectrum without CS 
interaction, iDessart fc Hillierl (120071 ) predict strong HV absorption in Hel 10830 A at a 
relatively late epoch; on day 50, the relative depth of Hel 10830 A is ~ 0.1 for an enhanced 
He abundance n(He)/n(He) = 0.2. This suggests a relative depth of ~ 0.05 for a solar 
He abundance. The absorption lies in the range of —12,000 to —17,000 km s~^, which 
corresponds to the blue wing of Hel 10830 A in SN 1999em on day 20 (Fig. [7]). The Hel 
10830 A line in SN 1999em on day 20 does not show evidence for the additional He absorption 
in the blue wing. Yet it is possible that in a more detailed model the time-dependent effects 
for He lines might be larger than in our simple model and would substantially contribute to 
the HV absorption of the He I 10830 A line at an age t > 20 d. However, for SN 1999em on 
day 20, we note the close correspondence between the w values deduced from the Ha and 
Hel 10830 A lines. 



Spectra of SN 1999em covering the He I 10830 A line are also reported by iHamuy et al.l 



( 200ll ) for three epochs: 8, 24, and 34 d. The first spectrum shows a strong Hel 10830 A 
line with a P Cygni proffie. This line, related to the early high temperature phase, rapidly 
disappears and is not seen in the next two spectra. Instead, HV absorption appears in this 
line. On day 24, the HV absorption of He I 10830 A is much the same as on day 20 in Fig. [TJ 
and it also appears on day 34. A steady HV absorption feature with time is expected in our 
model dlb)- On day 34, an additional redder absorption at ~ 10549 A appears separate from 
the HV absorption. The transformation of the He I 10830 A line from a single HV absorption 
into a double absorption structure is also observed between days 22 and 44 in another Type 
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IIP supernova, SN 1995V (iFassia et al.lll998l ). We believe that the low velocity absorption is 
the result of He e xcitation by gamrna-ray photon s leaking from the inner radioactive decays 
of ^^Ni and ^^Co (IFassia et al.lll998l : IChugailll987l ). To produce the required excitation of He 
about ~ 10~^ Mq of ^^Ni should be mixed out to a velocity ~ 4000 km s~^, which does not 



contradict the conclusion tha t the majority of the ^^Ni resides in deeper layers with velocity 
< 600 km s"^ Jutrobinlboorl ). 



4. ORIGIN OF THE NOTCH IN Ha 

We now address the issue of the notch that emerged in Ha in both SN 1999em and 
SN 2004dj at about day 60. We propose that the notch is produced by Ha absorption of the 
photospheric radiation in the CDS. We also consider an alternative mechanism: additional 
excitation of hydrogen in the unshocked ejecta by Lya photons emitted from the CDS. 



4.1. CDS vs. unshocked ejecta 

We first study the preferred explanation for the notch, Ha absorption in the CDS, and 
consider the clearly observed case of SN 2004dj. In modeling the Ha absorption produced by 
the CDS, the principal parameters are the column density of the cool gas, and its excitation 
temperature and turbulent velocity. The undisturbed or weakly disturbed CDS is charac- 
terized by a turbulent velocity of the order of the thermal velocity of the cool gas, ~ 10 — 20 
km s^^. However, with this turbulent velocity the equivalent width of the HV absorption 
turns out to be too small. A possible solution is to increase the turbulent velocity, as ex- 
pected by the action of the Rayleigh- Taylor (RT) instability of the decelerating thin shell. A 
significant deceleration of the model CDS of SN 2004dj, about 20% during the first 10 days 
(see low energy case in Fig. [T]), favors the development of the RT i nstability of the CDS 



followed by mixing of the coo l dense gas with hot forward shock gas (jChevalier fc Blondin 



19951 : iBlondin fc Ellisorul200ll ). The mixed cool gas thus can acquire a significant velocity 
dispersion and turbulent velocities up to ~ 0.1 of the shell velocity. Both components could 
contribute to the notch formation: absorption in the weakly disturbed CDS and in the mixed 
CDS material. We therefore consider two components: (1) a weakly perturbed spherical shell 
with a turbulent velocity of Ut ~ 20 km s^^ and (2) mixed CDS gas with a large turbulent 
velocity Ut ~ 400 — 500 km s^^. The mass fraction of the second component is denoted 
13. The corresponding spectral components of the HV absorption related to the CDS can 
be dubbed narrow notch and broad notch. These components should be distinguished from 
the 'ejecta component' - a broader shallow HV absorption formed in the unshocked ejecta 



- 13 - 



and addressed in the previous section. A cartoon (Fig. [8]) illustrates the origin the HV 
absorption produced by the ejecta and the CDS. 

Hydrogen excitation of the CDS gas is presumably maintained by the absorption of 
X-rays emitted from the reverse shock. The absorbed energy of X-rays is assumed to be 
homogeneously distributed through the CDS material, as expected for the modest shell 
column densities for our situation. The number density in the CDS is determined by pressure 
equilibrium P ^ p^v^ (where p„ is the pre-shock wind density) for the typical temperature 
of the cool gas of ~ 10^ K implied by the energy balance. On day 64, the CDS mass in 
our model for SN 2004dj is 3.2 x 10~^ Mq. The other input parameters provided by the 
interaction model are: CDS radius, = 5.2 x 10^^ cm, velocity Vg = 8350 km s~^. X-ray 
luminosity of the reverse shock Lx = 1.53 x 10^^ erg s~^, and electron temperature in the 
reverse shock kTx = 1.36 keV. The population of the second level of hydrogen is computed 
using a steady-state approximation (justified since Ue > 10^ cm~^), and assuming a two level 
plus continuum model for the hydrogen atom. Both non-thermal excitation and ionization 
by secondary electrons lead to the excitation of the second level, while depopulation proceeds 
by the two-photon transition, Lya escape and collisional deexcitation. 

The calculated profiles of HV absorption related to the CDS are overplotted on the 
model Ha profile with eje cta HV absorption and displayed in Fig. [9] along with the observed 



spectrum of SN 2004dj (IKorcakova et al.l |2005| ) on day 64. Shown are three cases: one 



without the CDS absorption and two cases with the CDS absorption that differ by the 
mass fraction of the broad CDS component (/? = 0.1 and 0.5). For the quoted values, the 
optical depths of the narrow components are 520 and 175, while these values for the broad 
component are 0.2 and 0.45, respectively. 

The absolute value of the radial velocity of the HV absorption produced by the CDS (fa) 
is slightly lower than the velocity of the CDS, Vg, owing to the sphericity effect. Assuming 
uniform brightness of the photosphere (Lambert's law) it is straightforward to show that the 
average absolute value of the radial velocity of the absorption is 

-a = ^^[l-(l-Xr^] , (7) 

where x = fp/fg is the ratio of the photospheric to the CDS velocity. According to this 
expression fa < Vg and, in the limit x <^ 1, fa = fg with high accuracy. The sphericity effects 
are responsible also for the additional broadening of the absorption over the radial velocity 
range {v^ — v^Y^'^ < \vr\ < Vs, which is ~ 460 km s~^ for the model of SN 2004dj on day 64. 

The alternative mechanism is that the notch could form as a result of an enhancement 
of hydrogen excitation in the unshocked SN ejecta due to the scattering of Lya photons 
emitted from the CDS that in turn is irradiated by X-rays from the reverse shock. A 
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resonance condition implies that Lya photons can scatter in the unshocked ejecta only in 
the narrow velocity range Vg < v < Vsn, where fgn is the boundary velocity of the unshocked 
SN ejecta. The expected flux of Lya from the CDS is 
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1 , -^X,abs 

2^ 47r R2 



where is the efficiency of Lya emission and Lx.abs is the X-ray luminosity absorbed by the 
CDS. The average excitation rate in the ejecta within the velocity width Av = Vsn — Vs is 
balanced by Lya escape; other depopulation mechanisms are negligible. The population of 
the second level (722) thus is approximately determined by the rate equation 



n2A2i/3u 
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where t is the SN age, huu is the Lya photon energy, and Pu is the Sobolev escape probability 
for Lya. For the model of SN 2004dj, the notch in Ha produced by Lya scattering in the 
ejecta is calculated on day 64 (Fig. [9]). Two cases, with efficiencies ip = 0.1 and 0.5, are 
plotted. One sees that even for the unacceptably high efficiency ip = 0.5 the notch produced 
by the Lya scattering is significantly weaker than that predicted by the absorption in the 
CDS. We conclude that the latter mechanism is the most likely explanation of the notch in 
Ha. 

We examined whether the absorption in the CDS might produce some effect in the He I 
10830 A line. For the above model, the optical depth in the He I 10830 A line turns out to 
be a factor of ~ 10^ lower than in Ha. Therefore, it is unlikely that the notch can be seen 
m He I 10830 A. The main reason is the efficient depopulation of the 2^S level in the CDS 
by Penning ionization of hydrogen. 

Given the relatively high Ha optical depth of the CDS, a notch at the same vel ocity 
could be seen in H/3 as well. We inspected spectra of SN 1999em (ILeonard et al.ll2002al ) and 
found the signature of a notch in H/3 a t the ri ght position —11,500 km s~^) between 
days 39 and 81. In fact, [Leonard et al.l (l2002al ) already discussed the presence of the HV 
absorption in Ha and H/3 in spectra of SN 1999em and concluded that these features cannot 
be caused b y metal line absorp tions. We also find a notch in H/3 in spectra of SN 2004dj 
obtained by IVinko et al.l (120061 ). Since the HV absorption components of Ha and H/3 have 
not been previousl y mentioned for this supernova, we show parts of the SN 2004dj spectra 
(jVinko et al.ll2006l ) retrieved from the SN archive SUSPECT (Fig. [TU]) . The spectra show 
conspicuous notches at about —8000 km s~^ both in Ha and H/3 between days 67 and 102. 
Below we will model HV absorption in both lines in more detail. 
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4.2. Rayleigh- Taylor mixing effects in HV absorption 

Here we present a model that can account for the strength and shape of the notch in 
Ha of SN 2004dj at the late photospheric epoch. Preliminary modeling suggests that the 
simple two-component model of the CDS (narrow and broad) introduced above is not able 
to describe the strength of the notch at late times (~ 100 d): the excitation produced by 
the X-ray absorption is insufficient. This suggests the presence of a more efficient excitation 
mechanism or an additional source of heat. Another problem is the relative depth of notches 
in the Ha and H/5 lines: they are comparable and shallow. This indicates that some veiling 
is present, possibly due to dumpiness. Both reasons force us to modify the two-component 
model of the CDS absorption. 

An important outcome of the RT mixing of the CDS material in the forward shock is the 
growth of the area of the contact surface {S) between cool and hot gas. The development of 
mixing may be thought as the progressive decrease of the minimum linear scale of fragments 
(A). This process is demonstrated in numerous experiments on the RT instability and in 
the case of SN/CS intera ction is clearly seen in 2D and 3D hydrodynamical computations 



(IBlondin fc Ellisonll200ll ). A fragmentation cascade of a fixed amount o f cool gas leads to 



the s caling 5* oc 1/A, while for a stationary mixing process S" oc A ( Sreenivasan et al. 



19891 ). The growth of S and, accordingly, of the surface-to- volume ratio for the mixed cool 
gas, should lead to a greater role of thermal conduction in the heating of the cool gas at late 
times. We speculate that this is the missing factor that could resolve the excitation problem. 

To describe the effect of thermal conduction, we introduce a third CDS component: 
an additional broad CDS component powered predominantly by thermal conduction. The 
mass fraction of the third component is yU -C 1, while the cumulative area of the contact 
surface for the third component, S = Cs^nR^, can be large, with an area ratio Cs ^ 1. We 
assume that the third component is homogeneously heated by thermal conduction, which is 
a sensible approximation for a small thickness of sheets of CDS gas, ~ 10^ cm, comparable 
to the mean free pass length of hot electrons and protons in the cool dense gas. The heat flux 
is presumably a small fraction of the saturated flux F^. = ?7c-fsat) 

we adopt 7]^ = 0.01 to allow 

qualitatively for the possible suppression of thermal conduction b y a magnetic fleld. Th e 



saturated flux is taken of the form Fgat = 0.4:{2kTe/7TmeY^'^nekTe ( Cowie fc McKed 119771 ). 



The forward shock is not equilibrated (Tg <^ Tj), so the ion (proton) conduction could be 
signiflcant. In our case, the ratio of proton-to-electron thermal flux is (Tj/Te)^''^(me/mp)^/^ ~ 
2 and this fact is taken into account. Ionization and thermal balance are solved for the dense 
hydrogen cooled by Lya and two-photon radiation. 

The RT mixing generally produces an inhomogeneous (lumpy) distribution of mixing 
regions. This is the result of the existence of a dominant angular scale of the RT instability. 
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about 10° — 15° fIChevalier fc Blondirull995l : lBlondin &: Ellisorul200ll ) and of the intermittence 



of the turbulent mixing. We describe the distribution of mixed regions as an ensemble of 
isolated lumps which occupy a fraction / ~ 0.1 — 1 of the volume; / = 1 corresponds to a 
homogeneous mixture of cool and hot phases. The fraction / should not be confused with the 
filling factor of the mixed cool dense gas, which is of the order of 10~^. Observationally, the 
lumpiness of mixing zones could be manifested as a veiling effect, i.e., a comparable moderate 
depth of absorption lines with different optical depth (viz.. Ha and H/3 as in our case). To 
facilitate the treatment of the radiation transfer in the lumpy medium we use an occultation 
optical depth, Tqc, i.e., the average number of lumps along the radius. This value is defined 
by / and a ratio of the total thickness of the mixing layer (forward post-shock layer AR) to 
the mean intercepted length of the lump d (for a sphere, d = |xradius): Tqc = f{AR/d). 
Generally, Tqc could be different for each of the three components. However, we assume 
the same value of Tqc for all three components. The effective optical depth in the line (r) 
then can be expressed through the average opt ical depth of homogeneous ly distributed CDS 



material and Toc- t = Toc[1 — exp(rav/Toc)] (IChugai fc Chevalierll2006l ). We thus consider 



a three component model of the cool dense gas which includes: (i) a narrow CDS component 
with turbulent velocity Ut ~ 20 km s~^, (ii) a broad CDS component with Ut ~ 400 — 500 
km s~^, both powered by X-rays, and (iii) a broad CDS component powered by the thermal 
conduction. 

The CDS notches in Ha and H/3 in SN 2004dj are calculated on the basis of the previous 
SN 2004dj model assuming ut = 20 km s^^ for narrow and Ut = 400 km s~^ for broad CDS 
components. The adopted fraction of the broad CDS component is (3 = 0.5 at the early 
epoch (t = 67 d), with a subsequent increase proportional to the age, while the fraction of 
the third component is /i = 0.1 (Table 2). The CDS absorption in Ha is superimposed on the 
model line profile with the ejecta HV absorption. For the H/3 line, the CDS HV absorption 
is superimposed on the background spectrum taken as a linear interpolation between the 
blue and red wings of HV absorption in the observed profile. The occultation optical depth 
Toe ~ 0.15 — 0.4 (Table 2) is chosen to satisfy the ratio Ha/H/3, while the area ratio Cs is 
determined from the strength of the absorption. Fig. ffTTl) shows that a three-component 
model with the area ratio Cs ~ 8 — 30 (Table 2) describes the HV absorption of both Ha 
and H/3 lines satisfactorily. On day 102 the ejecta HV absorption is somewhat stronger, 
which is responsible for the red-blue asymmetry in the model profile. However, this does not 
affect the intensity of the CDS absorption feature. The comparison to the two-component 
model and the behavior of Cs show that the role of thermal conduction increases with time, 
which is the expected result. Interestingly, the result depends rather weakly on fi. This is 
related to the canceling of effects of the variation of /i in the thermal excitation of hydrogen 
and the column density of the third component. We performed modeling of the CDS HV 
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absorption adopting the lower boundary velocity, Vc = 11,000 km s~^, and found that in 
this case deceleration of the CDS is too slow, which leads to a larger velocity of the CDS HV 
absorption on day 102 than is observed. The value Vc = 13, 000 km s"^ is about the least 
that provides the required deceleration rate of the CDS. 

Although the parameter values used for the third component are reasonable, the model- 
ing cannot be considered as proof that thermal conduction actually dominates the excitation 
of the CDS material at a late epoch - the model contains many assumptions. Yet the need 
for an additional source of excitation of hydrogen responsible for the HV absorption is an 
established result. Thermal conduction is a natural mechanism strongly suggested by the 
inevitable RT instability and mixing of the CDS matter with hot gas from the forward shock. 
On the other hand, the additional excitation might at least partially be related to enhanced 
soft X-ray radiation from hot dense inhomogeneities produced by the RT mixing and heated 
by thermal conduction up to ~ 



4.3. CDS absorption as a wind probe 



The identification of the features produced by Ha and H/5 absorption in the CDS pro- 
vides an efficient tool for the direct determination of the velocity of the SN/wind interface in 
SNe HP. With this value in hand we are able to obtain the wind parameter w using the thin 
shell model of the deceleration. Of course, the method assumes that the ejecta parameters 
M, E, and k are known. 

To facilitate the application of this method we use the self-similar solution. This de- 
scription is valid at the late photospheric epoch, when the initial conditions related to the 
structure of the outermost ejecta layers are already 'forgotten .' The self-simila r evolution of 
the C DS radius suggests R = At"^, with u = {k - 3)/{k - 2) Jchevahei]ll9823 jbl: iNadvozh 
19851 ). The factor A for the adopted density distribution (eq. [1]) is 



m 



A 



2Mu 



(fc-3) 



{k - 4:){k - 3)wC„ 



l/{k-2) 



(10) 



were Cm is defined by equation ([3]) and vq by equation ([2]). The parameters for SN 1999em 
are as before: M = 19 Mq, E = 1.3 x 10^^ erg, and k = 9. The radial velocity of HV 
absorption is calculated both for the Lambert photosphere (eq. [7]) and with limb darkening 
described by the first Ch andrasekhar approxini ation. The evolution of the photospheric 
velocity is set according to [Leonard et al.l (l2002al ). The observational estimat es of the radial 



velocity of HV absorption of Ha are obtained using the spectra reported by [Leonard et al. 
(l2002al ). The best fit for the wind density parameter is found tohe w = 0.99 ±0.09 (Fig. [T2 
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This estimate is in a good agreement with the value derived from the ejecta HV absorption 
of Ha {w = 1.00 ± 0.05) and the He I 10830 A hne {w = 1.15 ± 0.07). Uncertainties in the 
hmb darkening law and in the photospheric velocity do not affect the derived value of the 
wind density. 

In the case of SN 2004dj we use the same model as before: M = 19 Mq, E = 0.67 x 10^^ 
er g, and k = 9. The r adial velocities of HV absorption are measured from spectra obtained 



by lVinko et al.l (120061 ). The best fit wind parameter is w = 1.35 ±0.05 (Fig. [12]). This agrees 



with the value w = 1.27 ± 0.1 obtained for SN 2004dj from the ejecta HV absorption of Ha. 

The coincidence of w values found by diagnostics based on the HV absorption in ejecta 
and in the CDS is remarkable because the first probe is based on the model of X-ray emission 
by the reverse shock and the model of the ionization and excitation of ejecta, while the second 
is based only on the self-similar dynamics of the thin shell. 



5. DISCUSSION AND CONCLUSION 

Our goal was to study signs of CS interaction in the Ha and He I 10830 A absorption 
lines of SN IIP during the photospheric stage. We developed a model for the ionization and 
excitation of H and He in the unshocked ejecta taking into account time-dependent effects 
and the irradiation of ejecta by X-rays. The principal result is a prediction that for a typical 
RSG wind (w ~ 1), the effects of CS interaction in a normal SN IIP are significant and can 
be detectable at the photospheric stage in both lines. These effects consist of the emergence 
of: (i) HV absorption in the blue wing of the undisturbed Ha line (absorption 'shoulder') at 
t ~ 40 - 80 d, and (ii) HV absorption in He I 10830 A at an age t ~ 20 — 60 days at a radial 
velocity of ~ 10^ km s~^. We identify the HV absorption in SN 1999em in both Ha and He I 
10830 A and derive a wind density parameter w ^ 1. HV absorption in Ha is identified in 
SN 2004dj as well and the value of the wind density parameter is again found to be w ~ 1. 

The second important finding is that, in addition to the ejecta, absorption in the CDS 
formed at the SN/wind interface plays an important role in the HV feature, especially at the 
late photospheric phase. This component of the HV absorption is manifested as a notch in 
the HV absorption of Ha in SN 1999em and SN 2004dj after about day 60. A similar notch 
is seen in the H/3 line of these SNe. We developed a model that produces the notch in both 
Ha and H/3 at the different epochs of SN 2004dj which suggests a growing role for thermal 
conduction in the heating of the mixed CDS gas. The wind densities estimated from the 
observed velocities of HV absorption produced by the CDS in SN 1999em and SN 2004dj are 
consistent with the densities derived from the HV absorption produced by unshocked ejecta. 
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The ident i fication of HV absorption in Ha and H/3 in SN 1999em has been discussed by 
Leonard et aLl (j2002al ). Here, we identified HV absorption in Ha and H/5 in SN 2004dj at 
a velocity ~ 3000 km s~^ lower than in SN 1999em. The difference provides an additional 
argument against a possible metal line origin for these features. 



Leonard et al.l (j2002al ) claim that a similar component is also seen in the Nal doublet. 
Indeed, there is a weak dip with a depth of ~ 0.02 at the right velocity, ~ —11, 500 km s~^, 
if it is associated with the Nal 5890 A doublet. We calculated the optical depth of the Nal 
5890 A doublet for SN 1999em on day 60 in the model of absorption by the CDS and found 
that the narrow CDS component in this line has an optical depth of ~ 0.1, which means that 
after broadening by about 10 times the depth of absorption should become equal to ~ 0.01, 
comparable to the observed notch. However, we did not find a similar notch in the Na I line 
of SN 2004dj, or in any other SN IIP. At present, the reality of the HV absorption in the 
Nal 5890 A doublet remains uncertain. 

We examined other published spectra of SNe IIP at the photospheric epoch and foun d 
evidence for HV absorption in the Ha blue wi ng of SN 1985P (jChalabaev &: Cristianilll987l ). 
in both Ha and He I 1083 A of SN 1995V JPassia et al.l ll998h. in Ha of SN 2003gd. and 
SN 2OO60V Li et al.1 [20071 ). Available spectra of SN 1999gi do not cover the range 50 - 70 
days when the most pronounc ed CS interaction effect s should be observed. In noisy spectra 
on day 39 and 89 reported by [Leonard et al.l (l2002bl ). we do not find uiiequiy ocal evidence 
for HV absorption in Ha. We checked spectra of SN 2004et (ISahu et al.ll2006l ) and did not 
find convincing evidence for HV absorption in Ha or H/5. Early stage (20 — 40 d) spectra of 
SN 2004et show a shallow absorption in the Ha blue wing rapidly evolving towards the red. 
However the Sill 6347, 6371 A doublet may be responsible for this absorption. At a later 
epoch, a weak notch at 6282 A is present but unfortunately it coincides with the telluric 



absorption 6281.7 A line and the identification with HV absorption is doubtful. ISahu et al. 



(120061 ) consider the notch in Ha of SN 2004et to be real and remark that it might be related 
to CS interaction. 

If the HV absorption in SN 2004et is absent, the question arises why the dense wind sug- 
gested by a relatively h igh radio luminosity of this SN compared to SN 2004dj and SN 1999em 
(IChevalier et al.ll2006l ) is not revealed optically. The answer may be that the wind is so dense 
that the HV absorption is strong and gets merged with the main absorption as in the model 
case of w = 4 ( Fig. [21). This conj ecture is consistent with the strong H/? absorption observed 
in SN 2004et (ISahu et al.l l2006l ). A direct test of the proposed explanation might be the 
detection of strong He I 10830 A absorption in SN 2004et. Unfortunately, to our knowledge, 
observations o f the He I 10830 A region in SN 2004et are lacking. Interestingly, SN 2006my 
(ILi et al.l 120071 ). which does not show HV absorption, also has strong H/? like SN 2004et. If 



- 20 - 



the same mechanism (dense wind, w ^ 1) explains the absence of the HV absorption in this 
case, then SN 2006my should demonstrate high radio and X-ray luminosities for a SN IIP. 

We did not find signatures of HV absorption in Ha observations of th e low-luminosity 



SN 1 999br in spectra taken during the first 42 days after the explosion (jPastorello et al. 



20041 ). In the spectrum on day 42, the H/3 line is rather weak and the explanation that has 
been invoked for SN 2004et (strong CS interaction effect) is not applicable in this case. The 
likely reason for the absence of HV absorption in SN 1999br is a low wind density, w < 1. 
This conjecture has an interesting implicatio n. Since the mass loss rate i ncrea ses with the 
stellar mass according to the scaling law of iNieuwenhuijsen fc de Jagerl (Il990l ) , we expect 
that the main-sequence mass of pre-SN 1999br should be lower than that of SN 1999em. 

An interesting result, although model dependent, is that at the late photospheric stage 
of SN 2004dj [t ~ 100 d) thermal conduction may be the dominant heating mechanism of the 
mixed CDS gas in the forward shock. In this regard, we recall the Ha problem in SN 1979C: 
at an age of ~ 1 yr the H a luminosity was ~ 0.5 of the total reverse shock luminosity 



( IChevalier fc Franssonlll985l ) , which is beyond reasonable values of the Ha emission efficiency. 
A natural solution to the problem of the high Ha lumi nosity in SN 1 979C could be thermal 
conduction in the mixing layer of the forward shock (jChugail 119971 ). Another possibility, 
however, is that the mixing process could bring about enhanced X-ray emission from the hot 
dense mixed component, which may be another source of additional excitation of hydrogen 
in the CDS. This soft X -ray component ra ight explain the softening of the X-ray emission 
observed in SN 1999em JPoolev et al.ll2002h . 



Using Chandra X-ray data, iPooley et al.l (120021 ) conclude that the wind density param- 
eter for SN 1999em is in the range of w ~ 1 — 2 (or 1.5 — 3 for the distance 11.7 Mpc). Our 



present value 



w 



1) is within a factor of 1.5 of the lower limit. This should be considered 



as good agreement for independent estimates. For SN 2004dj the X-ray data are missing. 
The interpretation of radio data for SN 1999em and SN 2004dj in the model of wind free-free 
absorption indicates t hat the wind density in SN 2004dj is twice as low as in SN 1999em 
( IChevalier et al.ll2006l ). This is somewhat discrepant with the present results which indicate 
similar wind densities for these supernovae. The wind density of SN 1999em may be overes- 
timated compared to SN 2004dj from radio data because of the uncertain maximum in the 
8.47 GHz light curve o r a possible contribution of synchrotron self-absorption in SN 1999em 
(IChevaher etaPbood ). 



The wind density around SN 1999em and SN 2004dj suggests that the amount of the 
pre-SN material lost at the RSG phase (~ 10^ yr) is moderate in these cases. Assuming the 
wind velocity of a typical RSG, m = 10 km s~^, we find with w = 1 that at the RSG stag e 
the mass lost by the RSG wind is ~ 1 Mq. With the ejecta mass of 19 Mq (jUtrobinI 120071 ). 
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1.4 Mq enclosed in the neutron star, and < 1 Mq lost by the main-sequence wind, the initial 
mass of SN 1999em progenitor turns out to be about 22 — 23 Mq. 



The mass loss rate by massive stars is often estimated using the formula of lNieuwenhuijsen fc de Jager 



(ll990l ) in which M is defined through the stellar luminosity (L), mass (M), and radius (R). 
For the RSG stage, this formula can be modified using the relation betwe en the mass and lu- 
minos ity L/Lq = 14:{M / MqY'^'^'^ obtained from evolution calculations of iMeynet fc Maeder 
( I2OO3I ) in the range 15 — 25 for mode rate /no rotation. With this relation, the mass loss 
rate bv lNieuwenhuijsen fc de .Tageil f |l99nh becomes M = 2xlO-^\M/MQf-^^{R/RQ)^-^^ Mq 
yr~^. Since the luminosity at the RSG stage is determined by the mass of the He core which 
in turn is determined by the initial mass, we sh ould use the iii itial mass to estimate M. For 
SN 1999em with M = 22 Mq and R = 500 Rq dutrobinl [20071 ) this formula predicts a mass 
loss rate M ^ 8 x 10^^ Mq yr~^, i.e., eight times larger than our value. A similar result is 
expected for SN 2004dj. The disparity between our estimate of the pre-SN mass loss rate 
and the popular scaling law requires an explanation. The disparity cannot be attributed to 
a lower metallicity, because the metallicities of NGC 1637 (host galaxy for SN 1999em ) and 
NGC 2403 (host galaxy for SN 2004dj) are approximately solar jPilvugin et al1l2004l ). The 
difference can be reduced by a factor of 1.5 if we adopt a higher velocity for the pre-SN wind, 
M = 15 km s~^ instead of 10 km s~^. The choice of 15 kr n s~^ is suppo rted by the wind 
velocity of the massive RSG Betelgeuse, u = 14.3 km s~^ ( lMauronlll990l ). This, however, 
does not completely resol ve the problern . We admit the possibility that the hydrodynamical 
model of the light curve (lUtrobiru 120071 ) might significantly overestimate the mass of ejecta 
for some unknown reason. In that case, to obtain the wind w ^ 1 (for u = 15 km s~^) the 
initial mass of the presupernova must be M = 13.5 Mq with an ejecta mass of M 11 Mq 
instead of M ^ 19 Mq. With the SN energy reduced accordingly, the depth and position of 
the HV absorption remain the same (Fig. [2]c). 

To conclude, we have found three new probes of the wind density in SNe IIP that rely 
on spectroscopic observations during the photospheric epoch 20-100 d of (1) the ejecta HV 
absorption in Ha; (2) the ejecta HV absorption in He I 10830 A; and (3) the HV absorption 
(notch) produced by the CDS. The third method is more practical than the other two 
because it does not require modeling of the absorption lines; it uses velocity measurements 
and the self-similar model for the velocity fit. Of course, this method, as well as the other 
two, assumes that we know the ejecta parameters (M, E, k). However, the application of all 
three methods may be used to constrain ejecta parameters as well. The first two methods are 
sensitive to winds with a density parameter w > 0.5. However in Ha CS interaction effects 
saturate for large density w > 4. Therefore, only He I 10830 A remains a useful diagnostic 
tool for a high density wind with w > 2. For the third method the range of accessible 
wind density is less certain, because a CDS formed by the initial boundary shell could exist 
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in principle even in the case of a low density wind, although the hydrogen excitation and 
shell density depend on the wind density. Further understanding of the relation between the 
strength of deep HV absorption and the wind density requires 3D hydrodynamical modeling 
of the shock breakout phase in the presence of a wind and new spectroscopic observations 
in both Ha and He I 10830 A hues. 

We thank Daniela Korcakova for sending us spectra of SN 2004dj, Jozsef Vinko for 
informing us about SN 2004dj spectra in the archive SUSPECT, D. K. Sahu for sending us 
spectra of SN 2004et, and the referee for helpful comments on the manuscript. This research 
was supported in part by NSF grant AST-0307366. 
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Table 1: Parameters for the models of HV absorption 
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Table 2: Model parameters for the notch in SN 2004dj 
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Fig. 1. — Major properties of the interaction model for two values of the kinetic energy 
of the SN ejecta, 1.3 x 10^^ erg {left) and 0.7 x 10^^ erg. The top panel shows the X-ray 
luminosity of the forward [thick Ymes) and reverse {thin) shocks for unabsorbed {solid) and 
absorbed {dotte d) cases. In the lef t upper panel, Chandra X-ray luminosity measurements 



of SN 1999em (jPooley et al.l 120021 ) are plotted by crosses. The middle panel displays the 
electron temperature of the forward {thicklmes) and reverse {thin) shocks. The /ower panel 
shows the velocity of the thin shell {thick line) and boundary velocity of the unshocked SN 
ejecta {thin line). In both cases the X-ray luminosity is dominated by the reverse shock. 
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Radial velocity (10^ km s"^) 



Fig. 2. — The HV absorption in the Ha profile for different models (of. Table 1). Panel 
a shows the effect of wind density on day 50; panel b the evolution of the HV absorption 
in the standard model W on days 20, 50, and 80; panel c shows the neghgible effect of low 
ejecta density for the low mass case mW50; panel d shows the interaction effect in low energy 
SNe HP (cf. Table 1). 
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Fig. 3. — Schematic picture of the formation of Ha without and with CS interaction. In the 
absence of CS interaction (upper), the absorption component forms in the inner layers of 
ejecta (1) against the photosphere (Ph); the outer recombined ejecta (RE) do not contribute 
to the absorption hne profile (upper right). With CS interaction (lower), the double-shocked 
structure arises at the SN/wind interface with the forward shock (/), reverse shock (r) and 
contact surface where the cool dense shell occurs (c). The X-rays, primarily from the re- 
verse shock, ionize and excite the outer layer of ejecta (2) where the HV absorption forms, 
producing a depression in the blue wing of the undisturbed absorption (bottom right). 
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Fig. 4. — The same as Fig. 2 but for the He I 10830 A hne. 
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Fig. 5. — Modeling of the Ha profile with HV absorption in SN 1999em on days 41 and 54. 
The model {thick line) is overplotted on the observations (Elmhamdi et al. 2003). 
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Fig. 6. — Modeling of the Ha profile with HV absorption in SN 2004dj on days 55 and 64. 
The model {fMck line) is overplotted on the observations (Korcakova et al. 2005). On day 
64, a notch (n) appears against broader HV absorption (b). 
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Fig. 7. — High velocity absorption in the He I 10830 A line in SN 1999em on day 20. The 
model is overplotted on the observations (Elmhamdi et al. 2003). The thick hne is the 
model for w — 1.15 and the dashed line is the model with w — 1. Vertical bars show the rest 
positions of Paschen emission lines. 
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Fig. 8. — Visualization of Ha formation with a contribution from the CDS. Left: The 
cartoon of a SN depicts the photosphere {Ph.), the undisturbed atmosphere responsible 
for Hq; in the absence of CS interaction (i), the layer excited by X-rays from the reverse 
shock (2) responsible for the shallow ejecta HV absorption, the CDS {3n) responsible for 
the narrow notch, and the CDS gas mixed into the forward shock region (Sb) responsible 
for the broad notch. Right, upper to lower. The Ha profile without CS interaction (i), with 
the contribution of the ejecta HV absorption (1,2), with the contribution of a narrow CDS 
component {l,2,3n), and with the added broad CDS component {l,2,3n,3h). 
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Fig. 9. — Modeling of the HV absorption of Ha in SN 2004dj on day 64. The left panel 
shows, from top to bottom, the observed spectrum, a model without the CDS absorption, and 
models with the notch formed in the CDS assuming a two-component CDS with turbulent 
velocities of 20 and 500 km s~^ for different fractions of broad components (0.1 and 0.5). On 
the right panel, the two top lines are the same as on the left, while the next two are models 
of the notch formed in the unshocked ejecta by Lya absorption. The curves are labelled by 
the corresponding efficiency of Lya emission. 
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Fig. 10. — A sequence of SN 2004dj spectra (jVinko et al.l l2006l ) in the regions of the Ha 
and H/3 hnes for the same epochs, showing the HV absorption. Next to the Ha spectra are 
shown the corresponding ages; the ages of the H/3 spectra are the same. The notch is seen 
in both hnes at a radial velocity of about —8000 km s~^ starting on day 67. 
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Fig. 11. — Evolution of the HV absorption in the Ha and H/3 hnes of SN 2004 di in a model 



inclu ding absorption by CDS gas. The dotted lines are the observed spectra (IVinko et al. 



20061 ) ■ the thin line is the two-component CDS model composed of narrow and broad CDS 
components; the thick line is the three-component CDS model (Table 2) with the third 
component being a broad CDS component powered by thermal conduction. 
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Fig. 12.— Velocity of the HV absorption in Ha of SN 1999em and SN 2004dj. The crosses 
are the velocities of the HV absorption from observed spectra; the thick solid and dotted lines 
are the velocities of the HV absorption with and without limb darkening in the self-similar 
model; the thin solid line is the velocity of the CDS in the self-similar model. 



